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Spin interactions and switching in vertically tunnel-coupled quantum dots 
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KUngelbergstrasse 82, CH-4056 Basel, Switzerland 

We determine the spin exchange coupling J between two electrons located in two vertically tunnel- 
coupled quantum dots, and its variation when magnetic (B) and electric (E) fields (both in-plane 
and perpendicular) are applied. We predict a strong decrease of J as the in-plane B field is increased, 
mainly due to orbital compression. Combined with the Zeeman splitting, this leads to a singlet- 
triplet crossing, which can be observed as a pronounced jump in the magnetization at in-plane fields 
of a few Tesla, and perpendicular fields of the order of 10 Tesla for typical self-assembled dots. We 
use harmonic potentials to model the confining of electrons, and calculate the exchange J using 
the Heitler-London and Hund-MuUiken technique, including the long-range Coulomb interaction. 
With our results we provide experimental criteria for the distinction of singlet and triplet states and 
therefore for microscopic spin measurements. In the case where dots of different sizes are coupled, 
we present a simple method to switch on and off the spin coupling with exponential sensitivity using 
an in-plane electric field. Switching the spin coupling is essential for quantum computation using 
electronic spins as qubits. 
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I. INTRODUCTION 

Several methods to manipulate electronic spin in 
nanoscale semiconductor devices are being developed or 
are already available pi. Perhaps even more challenging 
is the proposal to use the electron spin in quantum dots 
as the basic information carrier (the qubit) in a quantum 
computer B- The recently measured long decoherence 
times in semiconductor heterostructures |3| and quantum 
dots 0] are encouraging for the further research of solid- 
state quantum computation. Quantum logic gates be- 
tween these qubits are effected by allowing the electrons 
to tunnel between two coupled quantum dots, thereby 
creating an effective spin-spin interaction. There is a 
large interest in quantum computation due to its poten- 
tial of solving some classically intractable problems, such 
as factoring Q, and speeding up the solution of other 
important problems, e.g. database search H. For the 
application of coupled quantum dots as a quantum gate, 
it is important that the coupling between the spins can 
be switched on and off via externally controlled parame- 
ters such as gate voltages and magnetic fields. In a recent 
publication M , we have calculated the spin interaction for 
two laterally coupled and identical semiconductor quan- 
tum dots defined in a two-dimensional electron system 
(2DES) as a function of these external parameters and 
have found that the interaction J can be switched on and 
off with exponential sensitivity by changing the voltage of 
a gate located in between the coupled dots or by apply- 
ing a homogeneous magnetic field perpendicular to the 
2DES. In this paper, we consider a different set-up con- 
sisting of two vertically coupled quantum dots with mag- 
netic as well as electric fields applied in the plane and 
perpendicular to the plane of the substrate (see Fig. |^). 
We also extend our previous analysis to coupled quantum 



dots of different sizes, which has important consequences 
for switching the spin interaction: When a small dot is 
coupled to a large one, the exchange coupling can be 
switched on and off with exponential sensitivity using an 
in-plane electric field Eu . 
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FIG. 1. (a) Sketch of the vertically coupled double quan- 
tum-dot system. The two dots may have different lateral di- 
ameters, aB+ and a_g~. We consider magnetic and electric 
fields applied either in-plane (By, E^) or perpendicularly {B±, 
E±). (b) The model potential for the vertical confinement is 
a double well, which is obtained by combining two harmonic 
wells at z = ±a. 

Semiconductor quantum dots are small engineered 
structures which can host single or few electrons in a 
three-dimensionally confined region. Various techniques 
for manufacturing quantum dots and methods for prob- 
ing their physical properties (such as electronic spectra 
and conductance) are known ||8|-|ll[]. In lithographically 
defined quantum dots, the confinement is obtained by 
electrical gating applied to a 2DES in a semiconductor 



heterostructure, e.g. in AlGaAs/GaAs. In vertical dots, 
a columnar mesa structure is produced by etching a semi- 
conductor heterostructure |12|. While laterally coupled 
quantum dots have been defined in 2DES by tunable 
electric gates [O-Om, vertically coupled dots have been 
manufactured either by etching a mesa structure out of 
a triple-barrier heterostructure and subsequently placing 
an electrical side-gate around it [^ or by using stacked 
double-layer self-assembled dots (SADs) |l^,|l^. In the 
mesa structure, the number of electrons per dot can be 
varied one by one starting from zero, whereas in SADs 
the average number of electrons per dot in a sample with 
many dots can be controlled, even one electron per dot 
is experimentally feasible [gOJ . 

Self-assembled quantum dots are manufactured in 
the so-called Stranski-Krastanov growth mode where a 
lattice-mismatched semiconducting material is epitaxi- 
ally grown on a substrate, e.g. InAs on GaAs |2^. Min- 
imization of the lattice mismatch strain occurs through 
the formation of small three-dimensional islands. Re- 
peating the fabrication procedure described above, a sec- 
ond layer of quantum dots can be formed on top of the 
first one. Since the strain field of a dot in the first layer 
acts as a nucleus for the growth of a dot in the second 
layer, the quantum dots in the two layers are strongly 
spatially correlated ||2^. Electrostatic coupling in verti- 
cal SADs has been investigated Il9], and it can be ex- 
pected that the production of tunnel-coupled SADs will 
be possible in the near future. 

In this paper, we concentrate on the magnetic prop- 
erties (including in-plane fields, B\\) of pairs of quantum 
dots in which two electrons are vertically coupled via 
quantum tunneling and are subject to the full Coulomb 
interaction. See Fig. |l| for a sketch of the system under 
study. Coupled quantum dots in the absence of quantum 
tunnelin g (p urely electrostatic interactions) were studied 
in Refs. (23|-|25|. Electronic spectra and charge densities 
for two electrons in a system of vertically tunnel- coupled 
quantum dots at zero magnetic field were calculated in 
Ref. [PGJ. Singlet-triplet crossings in the ground state of 
single |27[ and coupled dots with two [g^ to four [ p9|j30| 
electrons in vertically coupled dots in the presence of a 
magnetic field perpendicular to the growth direction {B± 
in Fig. Il|) have been predicted. 

In contrast to previous theoretical work on coupled 
dots [p3|-p0[ the investigation presented here both takes 
into account quantum tunneling and includes in-plane 
magnetic fields (i?|| in Fig. |^), leading to a much stronger 
suppression of the exchange energy than for B± (for 
very weakly confined dots, in-plane B fields can cause 
a singlet-triplet crossing, even in the absence of the Zee- 
man coupling). This result is in analogy to our earlier 
finding of a spin singlet-triplet crossing in laterally cou- 
pled identical dots as the perpendicular field is increased 
M. In addition to this, we investigate the influence of an 
electric field E± applied in growth direction on the low- 
energy electronic levels in the vertically coupled quan- 
tum dots. From the electronic spectrum, we derive the 



equilibrium magnetization as a function of both the mag- 
netic and the electric fields (magnetization measurements 
for many-electron double quantum dots are reported in 
Ref. |31|). As another important extension of earlier 
work, we consider a small dot which is tunnel- coupled 
to a large dot. We find that this system represents an 
ideal candidate for a quantum gate, since the exchange 
interaction J can be switched simply by applying an in- 
plane electric field E\\ (see Sec. ^. 

Our main interest is in the dynamics of the spins of the 
two electrons which are confined in the double dot. The 
spin dynamics can be described by an isotropic Heisen- 
berg interaction. 



Hs ~ J Si ■ S2 



(1) 



where the exchange energy J is the difference of the 
energies of the two-particle ground-state, a spin-singlet 
at zero magnetic field, and the lowest spin-triplet state. 
We shall calculate the exchange energy J(B, E, a) of two 
vertically coupled quantum dots containing one electron 
each as a function of electric and magnetic fields (E and 
B) and the inter-dot distance 2a. We show that an in- 
plane magnetic field has a much stronger infiuence on 
the spin coupling than a perpendicular magnetic field. 
Moreover, we will discuss the influence of the dot size 
on J, and investigate systems containing two dots of dif- 
ferent sizes. We will see that it is possible to suppress 
the spin-spin coupling exponentially by means of an in- 
plane magnetic field B\\ for large dots (weak confinement) 
or, alternatively, with an in-plane electric field Eu if one 
of the dots is larger than the other. Furthermore we 
will point out differences and similarities in the field- 
dependence of the tunnel-splitting t found in a quantum 
mechanically coupled double-dot system containing only 
a single electron and the exchange energy J, a quantity 
due to two-particle correlations. Performing these calcu- 
lations we make use of methods known from molecular 
physics (Heitler-London and Hund-MuUiken technique) 
thus exploiting the analogy between quantum dots and 
atoms. Note again that besides being interesting in its 
own right, a quantum-dot "hydrogen molecule", if ex- 
perimentally controllable, could be used as a fundamen- 
tal part of a solid-state quantum-computing device [0,0 , 
using the electronic spin as the qubit. 

In our discussion of the vertically coupled double-dot 
system we proceed as follows. In Section || we intro- 
duce a model for the description of a vertical double- 
dot structure. Subsequently (Sec. QI), we discuss ver- 
tically coupled quantum dots in perpendicular magnetic 
and electric fields. Sec. IV is devoted to the discussion of 



a double-dot structure in the presence of an in-plane mag- 
netic field. In Sec. M we present a simple switching mech- 
anism for the spin coupling involving an in-plane electric 
field. Finally, we discuss the implications of our result 
for two-spin and single-spin measurements in Sec. VI. 



II. MODEL 

The Hamiltonian which we use for the description of 
two vertically coupled quantum dots is 



H= Y^ h{r„pi) + C, 



/i(r,p)- — (p--A(r)) +ezE + Vi{r) + Vy{r), (2) 
C=^ r, 



K ri -r2 



where C is the Coulomb interaction and h the single- 
particle Hamiltonian. The dielectric constant k and the 
effective mass m are material parameters. The poten- 
tial Vi in h describes the lateral confinement, whereas Vy 
models the vertical double- well structure. For the lateral 
confinement we choose the parabolic potential 
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al+{x'^ +2/2), z >0, 
al_{x^ +y^), z <0, 



(3) 



where we have introduced the anisotropy parameters 
cko± determining the strength of the vertical relative to 
the lateral confinement. Note that for dots of differ- 
ent size (ao+ 7^ cto-) the model potential Eq. (0) is 
not continuous at z = 0. The lateral effective Bohr 
radii aB± = \/ fi / {muo zCio±) are a measure for the lat- 
eral extension of the electron wave function in the dots. 
In experiments with electrically gated quantum dots in 
a two-dimensional electron system (2DES), it has been 
shown that the electronic spectrum is well described by 
a simple harmonic oscillator pO|jll| ] . In the presence of a 
magnetic field i3j_ perpendicular to the 2DES, the one- 
particle problem has the Fock-Darwin states ||3^ as an ex- 
act solution. Furthermore, it has been shown experimen- 
tally iQ and theoretically [|2| that a two-dimensional 
harmonic confinement potential is a reasonable approxi- 
mation to the real confinement potential in a lens-shaped 
SAD. In describing the confinement V^ along the inter- 
dot axis, we have used a (locally harmonic) double well 
potential of the form (see Fig. [|b) 






(4) 



which, in the limit of large inter-dot distance a 3> cb, 
separates (for z « ±a) into two harmonic wells (one for 
each dot) of frequency uJz- Here a is half the distance 
between the centers of the dots and as = \/'hl(rnujz) 
is the vertical effective Bohr radius. For most vertically 
coupled dots, the vertical confinement is determined by 
the conduction band offset between different semiconduc- 
tor layers; therefore in principle a square-well potential 
would be a more accurate description of the real potential 
than the harmonic double well (note however, that the 
required conduction-band offsets are not always known 
exactly). There is no qualitative difference between the 



results presented below obtained with harmonic poten- 
tials and the corresponding results which we obtained 
using square- well potential s |34|| . 

It was shown in Refs. |7|,|35[| that the spin-orbit contri- 
bution (due to the confinement) i/so ~ {^l/'^meC^^ S • L 
with TTie being the bare electron mass can be neglected 
in the relevant cases, e.g. Hso/hujz ^ 10~^ for hujz ~ 
30 meV in GaAs. 

The Zeeman-splitting Hz ~ gi^-B X]i=i 2 -^ ' ^« ^^ ^'-'^ 
included in the two-particle Hamiltonian Eq. (0), since 
in the absence of spin-orbit coupling one can treat the 
orbital problem separately and include the Zeeman in- 
teraction later (which we will do when we study the low- 
energy spectra and the magnetization). Here, we have 
denoted the effective g-factor with g and the Bohr mag- 
neton with fiB- 



III. PERPENDICULAR MAGNETIC FIELD Bx 

We first study the vertically coupled double dot in a 
perpendicular magnetic field B — B± (cf. Fig. nl) which 
corresponds to the vector potential A(r) = B(— y, x, 0)/2 
in the symmetric gauge (for the time being, we set E — 
0). 

The confining potentials for the two electrons are given 
in Eqs. (Q) and (Q). As a starting point for our calcula- 
tions we consider the problem of an electron in a single 
quantum dot. The one-particle Hamiltonian by which we 
describe a single electron in the upper (lower) dot of the 
double-dot system is 



/»L(r)-^(p--A(r) 
Zm \ c 



+ - 



{al±{x^+y^) + {zTaf), 



(5) 



and has the ground-state Fock-Darwin 133] solution 



V±a(rj 



^ rnu!z \ 
\ nh J 



3/4 



-^^("±(^'+y') + (^=Fa)^)^ (g) 



corresponding to the ground-state energy e± — hujz{l + 
2a±)/2. In Eq. (^) we have introduced a±(B) ~ 

^al^+ujL{Bfliol = ^al^+B-^/Bl, with lol{B) = 

eB /2mc the Larmor frequency and Bq = 2mcu)z/e the 
magnetic field for which lOz = ojl- The parameters 
a±{B) describe the compression of the one-particle wave 
function perpendicular to the magnetic field. For find- 
ing the exchange energy J we make the Heitler-London 
ansatz, using the symmetric and antisymmetric two- 
particle wave-functions |*±) = (|12)±|21))/^2(1±5'2), 
where we use the one-particle orbitals (,c„(j(r) = (r|l) 
and (y9_|_a(r) = (r|2). Here \ij) = \i)\i) are two-particle 
product states and S = J d^rip*^^{r)ip^a{''^) = (2|1) de- 
notes the overlap of the right and left orbitals. A non- 
vanishing overlap S implies that the electrons can tun- 
nel between the dots. Using the two-particle orbitals 



|^±) we can calculate the singlet and triplet energy 
£s/t = (*±|^|*±) and therefore the exchange energy 
J = et — Es ■ We rewrite the Hamiltonian, adding and sub- 
tracting the potential of the single upper (lower) dot for 
electron 1(2) in H . The Hamiltonian then takes the form 
H = h'^^{ri) + h^^{r2) + W + C which is convenient, be- 
cause it contains the single particle Hamiltonians h'^^ and 
/i^a of which we know the exact solutions. The potential 
term is W^(ri,r2) = Wi{xi,yi,X2,y2) + W.u{zi, Z2), where 

Wi{xi,yi,X2,y2) = ^ Vi{xi,yi) 

1=1,2 

-^ H-i^l + vl) + (4+{A + vD] , (7) 



1=1,2 



W.{zi,Z2)= Y. V^{z,)-^[{zi+af + [z2-af]. 

(8) 



The formal expression for J is now 
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Evaluation of the matrix elements (12 |C + Vl^|12) and 
(12|C -t- M^|21) provides us with the resuh, 



J 



2^2 



rhuj. 



l-S'4 

c a-i- + a_ 



TT^l- {a++a^ - 1) 
1/9 9 \ a+ — a_ 

4 



V7Ie2M^(l-erf(dy2^)) 

arccos(a+ + a_ — 1) (10) 



+ 7 ("0+ ~ "0- 



;i-erf(d)) + -(l + d2) 



where erf(x) denotes the error function. We have in- 
troduced the dimensionless parameters d = a/ as for 
the inter-dot distance, and c = yJ-K /2{e^ / na-o) /fiuj^ 
for the Coulomb interaction. Note that a±, fi — 
2a+a_/ (a+ + a_), and the overlap 






(11) 



depend on the magnetic field B. The first term in the 
square brackets in Eq. ( pJ3| ) is an approximate evalua- 
tion of the direct Coulomb integral (12|C|12) for d > 0.7 
and for magnetic fields B < Bq |3^]. The second term 
in Eq. (10) is the (exact) exchange Coulomb integral 
(12|C|21)/S'2, while the last two terms stem from the 
potential integrals, which were also evaluated exactly. If 
the two dots have the same size, the expression for the 
exchange energy Eq. ( [lO| ) can be simplified considerably. 
We will first study the case of two dots of equal size, and 
later come back to the case of dots which differ in size. 

Setting ao+ = ao_ = uq in Eq. ( pi]| ) and using Eq. (11), 
we obtain 



J = 



fiojz 



sinh(2d2) 



zy^e^""'' (1 - erf (^^2^) ) (12) 



2a 



7r^l-(2a-l)2 



si2a-l) + -{l + d') 



where a ~ y^a^+B^/B^. As before, the first term in 
Eq. (O) is the direct Coulomb term, while the second 
term (appearing with a negative sign) is the exchange 
Coulomb term. Finally, the potential term in this case 
equals W = (3/4) (1 -I- d"^) and is due to the vertical con- 
finement only. For two dots of equal size neither the 
prefactor 25^/(1 — S*) nor the potential term depends 
on the magnetic field. Since the direct Coulomb term 
depends on B± only weakly, the field dependence of the 
exchange energy is mostly determined by the exchange 
Coulomb term. 

Note that for obtaining the large-field asymptotics 
{B > Bq) it would be necessary to include hybridized 
one-particle wave functions [[7| since in the magnetic field 
the level spacings between the one-particle ground states 
are shrinking and eventually become smaller than J, 
thus undermining the self-consistency of the one-orbital 
Heitler-London approximation. Increasing the inter-dot 
distance d (for fixed confinement hto), an exponential de- 
crease of the exchange energy J is predicted by Eqs. (nfl) 
and (|l2|). As mentioned, Eq. (|lC|) is an approximation 
and should not be used for small inter-dot distances 
d < 0.7. There are also some limitations to the choice of 
the anisotropy parameters ao± . If we consider a system 
with much stronger vertical than lateral confinement (e.g. 
ao± = 1/10), the exchange energy will become larger 
than the smallest excitation energy Ae = aa±hu;z in the 
single dot spectrum. In that case we have to improve our 
Heitler-London approach by including hybridized single- 
dot orbitals 0]. If, on the other hand, the two dots are 
different in size, double occupation of the larger dot is 
energetically favorable and a Hund-MuUiken approach 
should be employed. In the Hund-MuUiken approxima- 
tion, the Hilbert space for the spin singlet is enlarged by 
including two-particle states describing double occupa- 
tion of a quantum dot. Since only the singlet sector is 
enlarged it can be expected that we obtain a lower singlet 
energy eg than from the Heitler-London calculation (but 
the same triplet energy et) and therefore J = et — £s will 
be larger than the Heitler-London result, Eq. (|l^. 

We now apply the Hund-MuUiken approach to cal- 
culate the exchange energy of the double-dot system. 
We therefore introduce the orthonormalized one-particle 
wave functions $±a — {(p±a 
where g = {1 — yl 



5</'Ta)/Vl-2S'5-|- 



g' 



S^)/S. Using $±a, we generate four 
basis functions with respect to which we diagonalize the 
two-particle Hamiltonian H: The states with double oc- 
cupation, $^^(ri,r2) = <&±a(ri)<I'±a(i"2) and the states 
with single occupation ^!^(ri,r2) = [$+a(ri)$^a(r2) ± 
$_a(ri)$+a(r2)]/\/2. Calculating the matrix elements 
of the Hamiltonian H in this orthonormal basis we find 



H 



where 



f 2e + V+ -V2tn+ -V2iH- \ 

-V2tH+ 2e+ + U+ X 

-V2tH- X 2e- + U- 

V 2e + V- J 



e± = (*±a|/i(zTa)|$±a), e=-(e++e_), 



(13) 



(14) 



tH± = t - Z«± = -($±a|/l|<i>T.> - ^(ni^l^la), (15) 



V2^ +' 

y± = (*i|C|vI/^±), U± = (vl/d jc|*ij 
X={^iJC\^%). 



(16) 
(17) 
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FIG. 2. Left graph: Exchange energy J as a function 
of the magnetic field B applied vertically to the xy plane 
{B±, box symbols) and in-plane (S||, circle symbols), as cal- 
culated using the Hund-MuUiken method. Note that due to 
vertical orbital compression, the exchange coupling decreases 
much more strongly for an in-plane magnetic field. The pa- 
rameters for this plot correspond to a system of two equal 
GaAs dots, each 17 nm high and 24 nm in diameter (vertical 
confinement energy fiujz = 16meV and anisotropy parame- 
ter ao = 1/2). The dots are located at a center-to-center 
distance of 2a — 31 nm (d = 1.8). The single-orbital approxi- 
mation breaks down at about Bo ~ 9T, where it is expected 
that levels which are higher in the zero-field {B = 0) spectrum 
determine the exchange energy. Right graph: Single-particle 
tunneling amplitude t vs. magnetic field for the same system. 
Note that in contrast to the exchange coupling (a genuine 
two-particle quantity), t describes the tunneling of a single 
particle. Whereas J shows a weak dependence on the verti- 
cal magnetic field B±, we note that t{B±) (box symbols) is 
constant. 

The general form of the entries of the matrix Eq. (|l3) 
are given in Appendix |A|. The evaluation for perpendicu- 
lar magnetic fields B± can be found in Append ix O . We 
do not display the eigenvalues of the matrix Eq. (|l3|j here, 
since the expressions are lengthy. However, if the two 
dots have the same size (ao- = ao+), then the Hamil- 
tonian considerably simplifies since in- = te-i- = ^h, 
e+ = e_ = e and U+ = U- = U. In this case the eigen- 
values are e^i = 2e + C/h/2 + V^ ± V^h/4 + ^^h and 
Eso = 2e + [/h — 2X + V+ for the three singlets, and 
tt = 2t + V- for the triplet, where we have introduced 
the additional quantity Uh = U — V++X. The exchange 
energy is the difference between the lowest singlet and 
the triplet state. 



J = e,-e^^=V-^ + yUl + ml, (18) 

where we have used V — V- ~ V+. The singlet ener- 
gies es+ and Eso are separated from et and es_ by a gap 
of order Uu and are therefore negligible for the study of 
low-energy properties. If only short-range Coulomb in- 
teractions are considered (which is usually done in the 
standard Hubbard approach) the exchange energy J re- 
duces to -C//2 -f VC/2 + 16^2/2, where t and U denote 
the hopping matrix element and on-site repulsion which 
are not renormalized by interaction. We call the quan- 
tities in and C/h the extended hopping matrix element 
and extended on-site repulsion, respectively, since they 
are renormalized by long-range Coulomb interactions. If 
the Hubbard ratio tu /Uh is ^ 1, we are in the Hubbard 
limit, where J approximately takes the form (cf. Ref. |j) 



J = 



4t| 



V. 



(19) 



The first term in Eq. ( [191 ) has the form of the stan- 
dard Hubbard model result, whereas the second term V 
is due to the long-range Coulomb interactions and ac- 
counts for the difference in Coulomb energy between the 
singlet and triplet states ^j_. We have evaluated our 
result for a GaAs (to — 0.067TO,e, k — 13.1) system com- 
prising two equal dots with vertical confinement energy 
Tiu)z = 16meV (as = 17 nm) and horizontal confine- 
ment energy aQTiuJz = 8meV in a distance a — 31 nm 
{d = 1.8). The result is plotted in Fig. | (left graph, 
box-shaped symbols). The exchange energy J{B±) as ob- 
tained from the Hund-MuUiken method for two coupled 
InAs SADs (to = 0.08TOe |^, k = 14.6, Hlo^ = 50meV, 
ao+ = ao- = 1/4) is plotted in Fig. || (left graph, box 
symbols). Including the Zeeman splitting, we can now 
plot the low-energy spectrum as a function of the mag- 
netic field, see Fig. H (left). Note that the spectrum 
clearly differs from the single-electron spectrum in the 
double-dot (Fig. 0, right). 
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B(T) 
Exchange energy J (left 
gle-electron tunneling amplitude i (right graph) as a function 
of the applied magnetic field for two vertically coupled small 
(height 6nm, width 12 nm) InAs (m = O.OSme, n = 14.6) 
quantum dots (e.g. self-assembled dots) in a center-to-center 
distance of 9nm {d — 1.5). The box-shaped symbols corre- 
spond to the magnetic field B_i_ applied in z direction, the 
circle symbols to the field By in x direction. The plotted re- 
sults were obtained using the Hund-MuUiken method and are 
reliable up to a field Bo ~ 15 T where higher levels start to 
become important. 
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FIG. 4. Field dependence of the lowest four electronic 
levels for two vertically coupled InAs dots (parameters as 
in Fig. 0), including the Zeeman coupling with g- factor 
5inAs = —15. Left graphs (a,b): Spectrum for a two-electron 
system involving the Zeeman-split spin-triplet states (box, 
circle, and triangle symbols), and the spin-singlet (diamond 
symbols). The exchange energy J corresponds to the gap be- 
tween the singlet and the middle (m^ — 0, box-shaped sym- 
bols) triplet energy. Under the influence of an in-plane field 
B II a; (a), the ground state changes from singlet to triplet 
at about 9T, whereas in a perpendicular field B _L a; (b) the 
singlet-triplet crossing occurs at a higher field, about 12.5 T. 
Right graphs (c,d): Single-particle spectra, again plotted as 
a function of _B|| (c) and B± (d). Note that single-particle 
and two-particle spectra are clearly distinguishable. In par- 
ticular, there is no ground state crossing for a single electron. 
The B field dependence of the spectrum of the large GaAs 
dots (cf. Fig. bh is similar, with a much smaller Zeeman split- 
ting (gGaAs ~ —0.44). The plots are reliable up to a field 
-Bo ~ 15 T where higher levels start to become important. 

We now explain to what extent the Hund-MuUiken 
results (which we use for our quantitative evaluations 
of J) are more accurate than the results obtained from 
the Heitler-London method (which are more simple and 
which we used mostly for quahtative arguments). The 
Hund-MuUiken method improves on the Heitler-London 
method by taking into account double electron occu- 
pancy of the quantum dots. The Hubbard ratio tn /Ur 
can be considered a measure for the relative importance 
of double occupancy. Increasing the confinement TilUz at 
constant d (leading to potential wells that are deeper but 
closer together, since a = das = d^/hjmwl)^ we observe 
an increase in the discrepancy between Jhm and Jhl at 
zero magnetic field. Because the tunneling matrix ele- 
ment t is proportional to Tiujz and the on-site repulsion U 
is proportional to the Coulomb energy e^ / kub oc \/TujJz, 
the Hubbard-ratio tn /Uh increases as ^/huJ^ if the con- 
finement is increased at constant distance; thus double 
occupancy becomes more important, explaining the in- 
creasing difference between Jhm and Jhl- Both increas- 
ing the inter-dot distance 2a and the confinement hujz 
lead to a larger value of J = a/as and thus to a higher 
tunneling barrier. A strong decrease of the exchange 



energy J with increasing d is observed in both the re- 
sult calculated according to the Heitler-London and the 
Hund-MuUiken approach. 

We now turn to the dependence of the exchange en- 
ergy J on an electric field E± applied in parallel to the 
magnetic field, i.e. perpendicular to the xy plane. Using 
the Heitler-London approach we find the result 



J{B,Ej_) = J(B,0) + huJz 



2^2 



l-S'4 2 



E^ 
Eo 



o 1^ > (20) 



where Eg = muJl/eaB- The growth of J is thus propor- 
tional to the square of the electric field E±_ , if the field 
is not too large (see below). This result is supported 
by a Hund-MuUiken calculation, yielding the same field 
dependence at small electric fields, whereas if eE±a is 
larger than Uh, double occupancy must be taken into 
account. The electric field causes the exchange J at a 
constant magnetic field B to cross through zero from 
J{E = 0, B) < to J > 0. This effect is signalled by 
a change in the magnetization M, see Fig. g. 

In the presence of an electric field E±, the ground- 
state energy of an electron in the dot at z = ±a is 
e±(j;, B) = hcjz (1 + 2a±(B) - (J;/^o)' ± 2J^/^o) /2. 
The shift of the ground-state energies for the upper (e+) 
and lower (e_) dot due to an electric field can be used 
to align the ground-state energy levels of two dots of dif- 
ferent size (only for two dots of equal size, the energy 
levels are aligned at zero field). This is important be- 
cause level alignment is necessary for coherent tunneling 
and thus for the existence of the two-particle singlet and 
triplet states. The parameter Ea denotes the electric 
field at which the one-particle ground-states are aligned, 
e+{B,Ea) — e-.{B,Ea) (for dots of equal size, Ea = 0). 
Investigating the dependence of J on E±, one has to be 
aware of the fact that coherent tunneling is suppressed 
as the electric field is increased, since the single-particle 
levels are detuned (note, however, that the suppression is 
not exponential). This level detuning limits the range of 
application of Eq. (|2^) , which is only valid for small level 
misalignment, 2e{E± — Ea)a < J(0, 0), where J(0,0) is 
the exchange at zero field. Assuming gates at 20 nm be- 
low the lower and at 20 nm above the upper dot in the 
system discussed above (2a « 31 nm, Uujz = 16meV and 
ao = 1/2), we find that 2aE±e = J(0, 0) « 0.7meV 
at a gate voltage of about U w 1.6 mV. A further con- 
dition for the validity of Eq. (|2^) is J{Ej_) < hujzao^, 
(q^o- 5: Q;o+)- If this condition is not satisfied, we have to 
use hybridized single-particle orbitals. For the parame- 
ters mentioned above, we find J{E±) = huJzC^o- = 8 meV 
at a gate voltage U w 270 mV, therefore this condition is 
automatically fulfilled if 2eE±a < J(0, 0). The numbers 
used here are arbitrary but quite representative, as typi- 
cal exchange energies are on the order of a few meV and 
inter-dot distances usually range from a few nm to a few 
tens of nm. 

In the case where one of the coupled quantum dots is 
larger than the other, there is a peculiar non-monotonic 



behavior when a perpendicular field B± is applied at 
E = 0, see Fig. |5|. The wave-function compression due 
to the applied magnetic field has the eff'ect of decreasing 
the size difference of the two dots, thus making the over- 
lap Eq. (O) larger. This growth of the overlap saturates 
when the electron orbit in the larger dot has shrunk ap- 
proximately to the size of the orbital of the smaller dot, 
which happens at roughly Bq^ = 2mcuJzQ.o+/e (assuming 
that Q!o+ > ao-). 




FIG. 5. Exchange energy J as a function of the perpendic- 
ular magnetic field _B|| for two vertically coupled GaAs quan- 
tum dots of different size (both 25 nm high, the upper dot 
50 nm, the lower 100 nm in diameter, Bo+ ~ 2T). Here, J 
is obtained using the Heitler-London method, Eq. (lid). The 
non-monotonic behavior is due to the increase in the overlap, 
Eq. (|ll[), when the orbitals are magnetically compressed and 
therefore the size difference becomes smaller. 



IV. IN-PLANE MAGNETIC FIELD By 

In this section we consider two dots of equal size in 
a magnetic field -B|| which is applied along the x-axis, 
i.e. in-plane (see Fig. nl). Since the two dots have the 
same size, the lateral confining potential Eq. (^ reduces 
to V{x, y) — mujlaKx-^ + 2/^)/2, where the parameter ag 
describes the ratio between the lateral and the vertical 
confinement energy. The vertical double-dot structure 
is modeled using the potential Eq. (m. The single-dot 
Hamiltonian is given by Eq. (|5|) with the vector potential 
A(r) = B{0, —z, y)/2. The situation for an in-plane field 
is a bit more complicated than for a perpendicular field, 
because the planar and vertical motion do not separate. 
In order to find the ground-state wave function of the one- 
particle Hamiltonian /i+a, we have applied the variational 
method (cf. Appendix]^), with the result 



V± 



^« = (^j («o«/3) exp [±^^j (21 

'^ {aox^ + ay^ + (3{z T a)^) 



X exp 



Note that this is not the exact single-dot groundstate, 
except for spherical dots (ao = !)■ We have introduced 



the parameters a{B) — \/ 01^ -I- {B/BqY and /3(i?) = 
^/l -l- (B/Bq)^ describing the wave- function compression 
in y and z direction, respectively. The phase factor in- 
volving the magnetic length Is = \/hc/eB is due to 
the gauge transformation A±a = S(0, —{z =F cL),y)/2 -^ 
A = B{0,—z,y)/2. The one-particle ground-state en- 
ergy amounts to cq = ^^z (ao + a + (3) /2. From (p±a we 
construct a symmetric and an antisymmetric two-particle 
wave function $±, exactly as for B || z. Care has to be 
taken calculating the exchange energy J; Eq. (||) has to 
be modified, since ip±a is not an exact eigenstate of the 
Hamiltonian ft, j_^ (cf. Appendix y) . The correct expres- 
sion for the exchange energy in this case is 



J{B,d) ^ Jo{B,d)-huj, 



4S^ (3- 
l-S**" 



a 



Bo 



(22) 



where Jq denotes the expression from Eq. (0). The vari- 
ation of the exchange energy J as a function of the mag- 
netic field B is, through the prefactor 25^/(1 — 5^), de- 
termined by the overlap 



S{B, d) = exp 



PiB) 



1 



B 



a{B) \Bo 



, (23) 



depending exponentially on the in-plane field, while for 
a perpendicular field the overlap is independent of the 
field (for two dots of equal size), see Eq. (pTJ). We find 
that for weakly confined dots {hujz < lOmeV), there is 
a singlet-triplet crossing (J becoming negative), e.g. for 
huz = 7meV, Uq — 1/2, and 2a — 25 nm we find a 
singlet-triplet crossing at B « 6T. No singlet-triplet 
crossing is found for systems with vertical confinement 
fiWz ^ lOmeV. Generally, the decay of J becomes 
flatter as the confinement is increased. Improving on 
the Heitler-London result, we have again performed a 
molecular-orbital (Hund-MuUiken) calculation of the ex- 
change energy, which we plot in Fig. g (left graph, circle 
symbols). 
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FIG. 6. Magnetization M (in units of Bohr magnetons) 
as a function of the B field for vertically coupled large GaAs 
{g = —0.44) quantum dots (parameters as in Fig. 0) contain- 
ing two electrons (left graph) and a single electron only (right 
graph) at r = 100 mK. The box shaped symbols correspond 
to B±, the circles to By. The singlet-triplet crossing in the 
two-electron system (due to the Zeeman splitting and the de- 
crease of J) causes a jump in the magnetization around 5.5 T 
for _B||, but no such signature occurs for B±. 

It is crucial in experiments to distinguish between 
single- and two-electron effects in the double dot. A sin- 



gle electron in a double dot exhibits a level splitting of 
2t, where t denotes the single-particle tunneling matrix 
element, cf. Eq. (|l^), which has a B field dependence 
similar to the exchange coupling J. In order to allow 
a distinction between J and t, we have plotted t(B) in 
the right graph of Fig. |2[ Since the onc-particlc tun- 
neling matrix clement t is strictly positive, it is clearly 
distinguishable from the exchange energy J in systems 
with singlet-triplet crossing. Experimentally, the num- 
ber of electrons in the double-dot system can be tested 
via the field-dependent spectrum (Fig. ^) and magneti- 
zation (Figs. 
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FIG. 7. Magnetization Af (in units of Bohr magnetons) 
as a function of the B field for vertically coupled small InAs 
{g = —15) quantum dots (parameters as in Fig. M) con- 
taining two electrons (left graph) and a single electron only 
(right graph) at T = 4K. The box-shaped symbols corre- 
spond to B±, the circles to -B|| . The singlet-triplet crossing in 
the two-electron system causes a jump in the magnetization 
around 9T for B«, and at about 12.5 T for Bi. 
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FIG. 8. Magnetization M (in units of Bohr magnetons) as 
a function of the perpendicular electric field E_i_ for vertically 
coupled quantum dots containing two electrons at fixed mag- 
netic field. The box-shaped symbols correspond to B±, the 
circles to By . Starting at _E = with a triplet groundstate 
for _B|| (not so for -B±), the electric field eventually causes a 
change of the groundstate back to the singlet, which leads 
again to a jump in the magnetization for _B|| . The left graph 
corresponds to a GaAs double-dot (parameters as in Fig. Q) 
at T = 100 mK and B — 5T, whereas the right graph is 
for a smaller InAs double-dot (as in Fig. 0) at T = 4K and 
B = lOT. 



V. ELECTRICAL SWITCHING OF THE SPIN 
INTERACTION 

Coupled quantum dots can potentially be used as 
quantum gates for quantum computation [||J3], where 
the electronic spin on the dot plays the role of the 
qubit. Operating a coupled quantum dot as a quan- 



tum gate requires the ability to switch on and off the 
interaction between the electron spins on neighboring 
dots. We present here a simple method of achieving 
a high-sensitivity switch for vertically coupled dots by 
means of a horizontally applied electric field £'|| . The 
idea is to use a pair of quantum dots with different 
lateral sizes, e.g. a small dot on top of a large dot 
(ao+ > ao-, see Fig. |^). Note that only the radius in 
the xy plane has to be different, while we assume that 
the dots have the same height. Applying an in-plane elec- 
tric field E^^ in this case causes a shift of the single-dot 



orbitals by Ax± = eE» 
Eq = fiuj^/eaB-, see Fig. I 



2^,2 
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FIG. 9. Switching of the spin-spin coupling between 
dots of different size by means of an in-plane electric 
field _E|| [B = 0). The exchange coupling is switched 
"on" at _E = 0. When an in-plane electric field E\\ 
is applied, the larger of the two dots is shifted to the 
right by /S.X-, whereas the smaller dot is shifted by 
A2;+ < /S.X-, where Aa;± = E\JEoa^± and Eq = hcuz/eaB- 
Therefore, the mean distance between the electrons in 
the two dots grows as d' = \/<P + 'A^{E\\/Eo)'^, where 
A = (ao+ — ao-)/2cfo+ao_. The exchange coupling J, being 
exponentially sensitive to the inter-dot distance d' , decreases 
thus exponentially, J x S'^ ^ exp[-2yl^(iJ||/_Eo)^]. We have 
chosen huz = 7meV, d = 1, qo+ = 1/2 and oq^ — 1/4. For 
these parameters, we find Eo = hujz/eaB = Q.56mV/nm and 
A — (ao+ — ao_)/2Qo+Qfo_ = 6. The exchange coupling J 
decreases exponentially on the scale Eo/2A — 47mV//im for 
the electric field. 

It is clear that the electron in the larger dot moves 
further in the (reversed) direction of the electric field 
{Ax- > A.X+), since its confinement potential is weaker. 
As a result, the mean distance between the two electrons 
changes from 2d to 2d', where 



1 



d'= Jd^ + -{Ax^^Ax+y = Wrf2+^2 



l/a' 



0+ 



(24) 
)/2. Using Eq. (|ll|), we find 



with A = (l/c^o- 

that S oc exp(-d'2) oc exp[~'A'^{E\\/Eo)'^], i.e. the or- 



bital overlap decreases exponentially as a function of the 
applied electric field E\\ . Due to this high sensitivity, the 
electric field is an ideal "switch" for the exchange cou- 
pling J which is (asymptotically) proportional to S^ and 
thus decreases exponentially on the scale Eq/2A. Note 
that if the dots have exactly the same size, then A = 
and the effect vanishes. We can obtain an estimate of 
J as a function of -Ey by substituting d! from Eq. ( |24| ) 
into the Heitler-London result, Eq. (|lO|). A plot of J{E\\) 
obtained in this way is shown in Fig. |9| for a specific 
choice of GaAs dots. Note that this procedure is not ex- 
act, since it neglects the tilt of the orbitals with respect 
to their connecting line. Exponential switching is highly 
desirable for quantum computation, because in the "off" 
state of the switch, fluctuations in the external control 
parameter (e.g. the electric field ii^||) or charge fluctua- 
tions cause only exponentially small fluctuations in the 
coupling J . If this were not the case, the fluctuations in J 
would lead to uncontrolled coupling between qubits and 
therefore to multiple-qubit errors. Such correlated errors 
cannot be corrected by known error-correction schemes, 
which are designed for uncorrelated errors p^ ]. It seems 
likely that our proposed switching method can be realized 
experimentally, e.g. in vertical columnar GaAs quantum 
dots ||l^ with side gates controlling the lateral size and 
position of the dots, or in SADs where one can expect 
different dot sizes anyway. 



VI. SPIN MEASUREMENTS 

The magnetization (Figs. ^§|), measured as an ensem- 
ble average over many pairs of coupled quantum dots in 
thermal equilibrium, reveals whether the ground-state of 
the coupled-dot system is a spin singlet or triplet. On 
the one hand, such a magnetization could be detected by 
a SQUID or with cantilever-based [p8 39 1 magnetome- 
ters. This type of spin measurement was already sug- 



gested earlier for laterally coupled dots . The distinc- 
tion between spin singlet [S = 0) and triplet {S = 1) is 
also possible using optical methods: Measurement of the 
Faraday rotation (caused by the precession of the mag- 
netic moment around a magnetic field) using a pump- 
probe technique P,H] reveals if the two-electron system 
is in a singlet [S — 0) state with no Faraday rotation 
or in a triplet [S = \) state with finite Faraday rota- 
tion. Finally, it should also be possible to obtain spin 
information via optical (far-infrared) spectroscopy W%- 

We remark that if it is possible to measure the magne- 
tization of just one individual pair of coupled dots, then 
this is equivalent to measuring a microscopic two spin- 
1/2 system, i.e. 1/2 (g) 1/2 = ® 1. We have described 
elsewhere how such individual singlet and triplet states 
in a double dot can be detected (through their charge) 
in transport measurements via Aharonov-Bohm oscilla- 
tions in the cotunneling current and/or current correla- 
tions |^-||. 



It is interesting to note that above scheme allows one 
to measure even a single spin 1/2, provided that, in addi- 
tion, one can perform one two-qubit gate operation (cor- 
responding to switching on the coupling J for some finite 
time) and a subsequent singic-qubit gate by means of ap- 
plying a local Zeeman interaction to one of the qubits. 
[Such local Zeeman interactions can be generated e.g. by 
using local magnetic fields or by inhomogeneous g- factors 
pO| .] Explicitly, such a single-spin measurement of the 
electron is performed as follows. We are given an arbi- 
trary spin 1/2 state \a) in quantum dot 1. For simplicity, 
we assume that |a) is one of the basis states, |a) = 1 1) 
or \a) = II); the generalization to a superposition of the 
basis states is straightforward. The spin in quantum dot 
2 is prepared in the state | "f). The interaction J between 
the spins in Eq. (nl) is then switched on for a time Ts, 
such that Jp" J{t)dt = it /A. By doing this, a 'square- 
root-of-swap' gate [0,pHl is performed for the two spins 
(qubits). In the case \a) = \]), nothing happens, i.e. the 
spins remain in the state 1 1|), whereas if \a) = | i), then 
we obtain the entangled state (| J,t) +«l Ti))/^/2i (up to a 
phase factor which can be ignored). Finally, we apply a 
local Zeeman term, g^gBS^, acting parallel to the z-axis 
at quantum dot 1 during the time interval tb, such that 
Jo^ {9f^BB){t)dt = 7r/2. The resulting state is (again up 
to unimportant phase factors) the triplet state | IT) in 
the case where \a) = 1 1)1 whereas we obtain the singlet 
state (I ti> - I iT>)/\/2 in the case |a) = U). In other 
words, such a procedure maps the triplet | ff) into itself 
and the state | 11) into the singlet (similarly, the same 
gate operations map \ll) into itself, while | "f j) is mapped 
into the triplet (| fi) -I- | it))/V2j again up to phase fac- 
tors). Finally, measuring the total magnetic moment of 
the double dot system then reveals which of the two spin 
states in dot 1, | f) or 1 1), was realized initially. 



VII. DISCUSSION 

In summary, we have calculated the spin exchange in- 
teraction J(B,E) for electrons confined in a pair of ver- 
tically coupled quantum dots, and have compared the 
two-electron spectra (with level splitting J) to the single- 
electron spectra (with level splitting 2i) . Comparing the 
one- and two-electron spectra enables us to distinguish 
one-electron filling from two-electron filling of the dou- 
ble dot in an experiment. For two-electron filling in the 
presence of a magnetic field, a ground-state crossing from 
singlet to triplet occurs at fields of about 5 to 10 T, de- 
pending on the strength of the confinement, the coupling, 
and the effective g-factor. The crossing can be reversed 
by applying a perpendicular electric field. 

As a model for the electron confinement in a quantum 
dot we have chosen harmonic potentials. However, in 
some situations (especially self-assembled quantum dots) 
it is more accurate to use square- well confinement poten- 
tials in order to model the band-gap offset between differ- 



ent materials. We have also performed calculations using 
square-well potentials, which confirm the qualitative be- 
havior of the results obtained using harmonic potentials. 
The results from using the square-well model potentials 
cannot be written in simple algebraic expressions, and 
are given elsewhere [Q . 

Furthermore, we have analyzed the possibilities of 
switching the spin-spin interaction J using external pa- 
rameters. We find that in-plane magnetic fields B\\ (per- 
pendicular to the inter-dot axis) are better suited for 
tuning the exchange coupling in a vertical double-dot 
structure than a field B^ (applied along the inter-dot 
axis). Moreover, we have confirmed that the dependence 
of the exchange energy on a magnetic field is stronger for 
weakly confined dots than for structures with strong con- 
finement. An even more efficient switching mechanism is 
found when a small quantum dot is coupled to a large 
dot: In this case, the coupling J depends exponentially 
on the in-plane electric field E\\ , and thus provides an 
ideal external parameter for switching on and off the spin 
coupling with exponential sensitivity. The experimental 
confirmation of the electrical switching effect would be 
an important step towards solid-state quantum compu- 
tation with quantum dots. 

Another (very demanding) key experiment for quan- 
tum computation in quantum dots is the measurement 
of single electron spins. We have presented here a 
theoretical scheme for a single-spin measurement using 
coupled quantum dots. Obviously this scheme already 
requires some controlled interaction between the spins 
(qubits) and therefore the successful implementation of 
some switching mechanisms. 



X = N^ [(1 + g^)G\S^ + .g2(G+ + G") + 2g^S^G2 

-f252G3-2g(l + g2)5(G+ + G4)], (A4) 
w± = N^ [~gGf - g'^Gf - .g(l + g^){2S^G'i + G3) 

+5(f + 3g2)G± + 5V(l + ff')Gj], (A5) 



with TV = 1/ VI - 2Sg + g2 and 5 = (1 - Vl-S^)/S. 
We have introduced the overlap integrals 



Gf = {ip±a'P±a\C\(p±aV±a), 
G? = S~'^{ip±a(p±a\C\f^a'PTa), 
G2 = S^^{ip±a(p^a\C\ipzfzaip±a), 
G3 = {'P±aVTa\C\V'±aVTa) , 

Gf = S^'^ {ip±aip±a\C\ip±aip^a) ■ 



(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 



Note that the expressions for G?, G2, and G3 are invari- 
ant under exchange of ipa and (p-a- In the case where 
the two single-particle Hamiltonians coincide (implying 
that the dots have the same size), we find Gf = G^ 
{— Gj, since C depends only on the relative coordinate) 



4 , and the expressions in Eqs. (A1)-(A5) 



and GJ = G 

for the matrix elements can be simplified accordingly. 
This simplification leads to the same form of the Hund- 
MuUiken matrix elements which we have calculated for 
laterally coupled dots jj]. If it is possible to choose the 
orbitals ip±a to be real, e.g. if the magnetic field is in z 
direction, then G^ = G2, leading to a further simplifica- 
tion of the matrix elements Eqs. (Al)-(A5). 



APPENDIX B: HUND-MULLIKEN MATRIX 

ELEMENTS, B ±X,Y 
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APPENDIX A: HUND-MULLIKEN MATRIX 
ELEMENTS 

Here we list the explicit expressions for the matrix el- 
ements defined in Eqs. (13)-(p^) for two dots with ar- 
bitrary (and possibly different)single-electron Hamilto- 
nians h±a and (non-orthogonal) single-electron orbitals 
ip±a centered &t z = ±a. The matrix elements are 



V+ = N^ [2g^(G+ + G]:)+ Ag^S'^G^ + Ag'^G; 

+ (l + g2)2G3-6g2(G+ + G: 
V^^N\l~g''f[G3^S''G2], 
U± = iV" [Gf + g^Gf + 2g^S^Gl + 2g^S^G2 

+2g^S^G3 - 4gS{Gt - g^Gf)] , (A3) 



(Al) 
(A2) 



If the single-electron Hamiltonian is given by Eq. (H) 
with a perpendicular field B J- x,y then we can further 
evaluate the integrals Eqs. (A6) -(|A10| ) and the single- 
particle matrix elements in Eqs. (|13|)~(|I7|) as a function 
of the dimensionless inter-dot distance d = o/ob and the 

B'/Bi. 



magnetic compression factors a±{B) — y uq_j_ 
The single-particle matrix elements are given by 



hco. 



e± 



1 



16d2 1 - 5*2 V 5 



a± 



ga^ 



±- 



1 a, 



0+ '-'o- 



4 a+a- 



ga± I (1 — erf(rf)) 
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^' ^- (I + <f) 



1 - S'2 V 4 



{a± +q;=p) 



(Bl) 



TlLJz S 

"1-5*2 



ic4+_(4_ 

4 «+«_ 



(l-erf(d))(a+-a-) 



+ \{l + d') 



(B2) 



where we have used S — \^^/a+aI / {a+ + a-)] exp(— fi2). 
The (two-particle) Coulomb matrix elements can be 



10 



expressed as in Eqs. (A1)-([A5|), where the integrals 
Eqs. ([A^)- (|A10| ) take the form 



Gf == Hlu, 



2c 



a± 



^ Vl-(2a±-l)2 



arccos(2a± — 1), 



t- c (a+ + a_) arccos(a+ + a_ — 1) 
G2 = Ct]^ = ncuz — - 



v/1 - (a+ +a- 
G3 = ^Jcu^cV^e^M'^' (1 - erf(dv^)) , 

G4 



1)^ 



(B3) 

(B4) 
(B5) 



;ic^,cW ^",^^"+^"'^ e^^'^'[l-erf(dVMl)], (B6) 



So-i 



xe K/5+"o)'-=-ia!/= cos ( 2yd^ ] , (C5) 



G4 = G J = G4 = hiVz ■^— 5- \/aao/3 



OC pQO 



dy dzKo(^{y' + z')) 



00 J —00 



xe 



-\(2a-ao)y^-^(3(z-d)^ + \aoz 



B 



cos (yd—). (C6) 



Here Kg denotes the zeroth order Macdonald function 
and Iq is the zeroth order modified Bessel function. The 
quantities a, j3 and S have been defined earUcr. 



where we have introduced /i — 2a+a_/(a+ + q; _) a nd 
^± = {a\ + a+a_) / (3a± + a^). Note that Eq. (B5) is 
an approximation which deviates from the exact result 
by less than 12% in the range d > 0.7 and /it < 1 as we 
have checked by numerical evaluation of the integrals. 



APPENDIX C: HUND-MULLIKEN MATRIX 

ELEMENTS, B\\X 

The Hund-Mulliken calculation for a system of two 
equal dots with a magnetic field applied in x-direction 
(Sec. [V) is formally identical to the one with a field in 
z-direction presented in Sec. Ill, For equal dots we set 

a, and e+ = e_ = e. The 



ao+ = ao- = ao, a+ = a. 
one-particlc matrix elements are then 






hiOz 



ao + a + f3 + 



52 



16d2/32 
1-52 



1-524 V/3 



- - + d^ 



2d 



2 1-52 



3/1 



2d 



,2 /? - a ^2 
a B^ 

,13-aB^' 

^rBi_ 



(CI) 
(C2) 



Since we consider two equal dots, the matrix elements 
of the Coulomb Hamiltonian are formally equal to the 
matrix elements given in Ref. M, where Fi has to be 
replaced by d defined by 



Gi = G+ = G^ = G? = huj.-^/aaoP / drr 

Ti" Jo 



c 



G2 — huJz — \/aaof] dr dz 

" Jo J-oc 



\/r^ + z2 



X lo ( " ~ "" r2 1 e-K"+"o)r=~i/3(2+2d)=^ 



(C3) 



(C4) 



c 



G3 = ?ic^.3v/^^e^(^/^«) / 
dy 



X / dr 
10 



\Jr'^ + y- 



..l,\^^^r' 



APPENDIX D: HEITLER-LONDON 
CALCULATIONS, B \\ X 

In the following we evaluate the exchange energy J 
for two coupled quantum dots in a magnetic field ap- 
plied perpendicularly to the inter-dot axis (B||x) using 
the Heitler-London approach. We first study the one- 
particle problem for an anisotropic quantum dot with a 
magnetic field applied perpendicularly to the symmetry 
axis of the dot, 



/i°(r 



(Dl) 



where ao is the ellipticity and A(r) — B{0,—z,y)/2. 
We can separate /i°(r) = h'^ix) + hy^{y,z) into a B- 
independent harmonic oscillator 



^2 n2 



,r,, , h d muj; n 2 



(D2) 



and a _B-dependent part 



mzLU / 2 2 , a2 2\ 



hl{y^z)^pl+pl-u,L^+'^{a^y^+P^z^), 



(D3) 



with a = ^af, + {u jLluJz Y = ^af^ + {B/BoY, and 
/3 = yJl + iujL/uj,Y = y/T+jBjBof- We have not 



solved Eq. (D3) exactly; instead we have used a varia- 
tional approach, minimizing the single-particle energy 



eo 



(^1^) 



(D4) 



as a function of two variational parameters, in order to 
find a good approximate ground state wave function. A 
reasonable trial wave function ip should reproduce the 
anisotropy between y and z in the Hamiltonian. This 
requirement is fulfilled e.g. by a Gaussian 



V'i(7i,72,y,z) =Mt 



— llV —722 



(D5) 
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or by mixing Fock-Darwin states V'o,; with angular mo- 
mentum / = 0, 2, —2 and radial quantum number n — 0, 



tp2{S2,S-2,y,z) =Af 



i^o,oiy,z) 



E 

l=±2 



SiipoAy^^) 



(D6) 



where 5-2, 62, and 71,72 are variational parameters 
and Af, M are normalization constants. Calculating 
60(71,72) and eo(5-2,'^2), and subsequently minimizing 
with respect to the variational parameters, we find that 
ilii{mijJza/ {2fi),muizl3 / {2h),y, z), with the normalization 
constant M — {mujz/'n'hy^'^{apy^^ is the best approxi- 
mate ground-state wave- function in our variational space. 
We have also shown that including the Fock-Darwin 
states with angular momentum quantum numbers ^ = ±1 
in ■02 does not lead to a lower minimum of the energy 
(V'2|^yzl'02)/(V'2|''/'2)- The full one-particle wave function 
is then given by 

(D7) 

Shifting the single-particle orbitals to (0,0, ±a) in the 
presence of a magnetic field we obtain ^±a{x,y,z) = 
exp(±iya/2^^)(^(a;, y,z^ a). The phase factor involving 
the magnetic length Ib = ^Uc/eB is due to the gauge 
transformation A±a ~ B{0, — {z^ a) ,y)/2 -^ A — 
-6(0, — 2, y)/2. Having found an approximative solution 
for the one-particle problem in a dot centered at z = +a 
or z = —a, we show that the exchange energy is given by 
Eq. (P^ ) for a system with two dots of equal size, where Jo 
denotes the result from Eq. (H). In the derivation of the 
formal expression for the exchange energy Jo{B, d) given 
in Eq. (^) , we have used that Lp±a was an exact eigenstate 
of ft,[j_^, and therefore ((p=pQ|/i±„|v9±a) = S'(<p±a|/i±a|(/'±a), 
where S — {ipa\yy-a) denotes the overlap of the shifted 
orbitals. The approximative solution Eq. (D7) for an 
anisotropic dot in the presence of an in-plane magnetic 
field is not an exact eigenstate of /i". Using the corrected 
off-diagonal matrix element 



{V>Ta\h±a\V>±a) ^ S 



— ^ ao + a + /3 +cfV^ 



(D8) 



the result for the exchange energy Eq. (22) can easily be 
derived. 



* Electronic address; guido.burkard@unibas.ch 
^ Electronic address: daniel.loss@unibas.ch 
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